Interface dependence of the Josephson-current fluctuations in short mescoscopic 
superconductor / normal-conductor / superconductor junctions 
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We discuss the dependence of the Josephson-current correlations in mesoscopic 
superconductor/normal-conductor/superconductor (SNS) devices on the transparency of the 
superconductor/normal-conductor interfaces. Focusing on short junctions we apply the super- 
symmetry method to construct an effective field theory for mesoscopic SNS devices which is 
evaluated in the limit of highly and weakly transparent interfaces. We show that the two-point 
Josephson-current correlator differs by a universal factor of 2 in these two cases. 

PACS numbers: 74.45.+c,74.50.+r,73.23.-b 
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I. INTRODUCTION 

One of the most promising directions in the field of 
superconductivity in mesoscopic structures is provided 
by recent technological advances in hybrid supercon- 
ductor/semiconductor technology. As an example we 
want to mention the very recently observed quantiza- 
tion of the critical current in superconducting quan- 
tum point contacts^. An important aspect of this de- 
velopment is that for hybrid semiconductor structures 
both the system geometry and the chemical potential 
can be varied flexibly, which greatly facilitates the ob- 
servation of mesoscopic fluctuations. In contrast, for 
ordinary metal heterostructures only external magnetic 
fields effectively qualify as averaging parameters, imply- 
ing that the dependence of fluctuation phenomena on the 
phase of the superconducting condensate is difficult to 
observe. Indeed, fabrication of such hybrid devices has 
advanced to the stage that fluctuation phenomena may 
be measured^iiiSi&l in a wide range of geometries and 
conditions. 

The observable behavior of a two-dimensional electron 
gas (N), as realized in a semiconductor herterostructure, 
in contact with a superconductor (S) is highly sensitive to 
the quality of the SN contact. In this work we discuss the 
SN interface dependence of the Josephson-current J{4>) 
through an electron gas, sandwiched between two super- 
conductor terminals whose order parameters exhibit a 
phase difference <fi. A measurement of the fluctuation 
behavior var J{4>) of this quantity has been attempted 
already by Takayanagi et al£. As pointed out above, the 
use of a semiconductor heterostructure allows the use 
of an external gate voltage to tune fluctuations in the 
supercurrent, while keeping the phase difference fixecf . 
For typical experiments^ transmission tends to be weak 
due to the presence of the Schottky barrier. However, as 
noted in Ref. 0] , a theory of supercurrent fluctuations in 
the weak transmission regime has so far been missing. In 
this work, we fill this gap and compare the behavior of 
the supercurrent for weakly and highly transparent in- 
terfaces. (We remark that an experimental realization of 
a crossover in the transition strength has already been 



attempted in Ref.Q.) Specifically, we will compute the 
mesoscopic fluctuations of the current, var J(<fi) and, more 
generally, its correlation profile 



Kfafa) = (J(fa)J(<h)) - (J(^))(J(h))- 



(1) 



This latter function provides a wealth of information re- 
lating to the design of the system, and it is clear that its 
experimental determination - along the lines of Ref. Q - 
would provide a sensitive test of the theory of mesoscopic 
fluctuations in SN systems. However, before proceeding, 
let us briefly put our present analysis into the context of 
previous studies of the problem. 

A pioneering work in the theory of supercurrent fluctu- 
ations is due to Altshuler and Spivak^, who calculated the 
supercurrent correlations by means of a diagrammatic 
perturbation theory. Their results assumed the limit of a 
long junction, so that E c A, where E c is the Thouless 
energy and A the superconducting order parameter. In 
this limit, the supercurrent fluctuations are not univer- 
sal: instead, the supercurrent variance is proportional to 
{eE c /K)\ 

The applicability of the diagrammatic theory is lim- 
ited, however, to conditions under which the fluctuations 
of the supercurrent exceed its average. While such a sit- 
uation may be reached by application of a magnetic field, 
or in the presence of a glassy structure in the N region, 
it clearly is rather specialized. Under more usual condi- 
tions, the essentially non-perturbative nature of the prox- 
imity effect leads to a proliferation of diagrams whose 
summation is practically impossible^. 

A powerful alternative to the diagrammatic approach 
is provided by a multiple scattering theor y 11 ! 12 . Here 
the property in question is related to the scattering ma- 
trix of the N region, whose statistical properties are 
knownii. The variance of the supercurrent through the 
N region at fixed phase difference has been calculated by 
Beenakke r 11 ! 12 using the scattering approach. His results 
apply for a short junction A -C E c with highly transpar- 
ent SN interfaces. In this case the fluctuations become 
universal: the supercurrent variance is proportional to 
(eA/h) 2 . 

In this work, we apply the supersymmetry method 13 
to construct an alternative approach to the problem. A 
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general feature of our formulation is its close alliance to 
the quasiclassical approac h 10 ' 14 i 15 ' 16 , itself a traditional 
means of describing the mean properties of SN systems. 
In particular, using the fact that the stationary phase 
configurations of the field theory of dirtyii SN systems 
are determined by the quasiclassical Usadel equations 10 , 
one can benefit from the huge body of expertise on (en- 
semble averaged) Josephson-currents in mesoscopic SN 
devices. The evaluation of fluctuations around the Us- 
adel mean field then leads to results for the mesoscopic 
fluctuations of the current. 

To be specific, we apply our approach to the analysis 
of short junctions with moderately weak coupling of the 
superconductor condensate to the normal-metal "quan- 
tum dot" , which allows us to use random-matrix theory 
(RMT) to model the normal metal. That is, we assume 
the hierarchy of energy scales A <C E g <C E c , where 
Eg = dg is the inverse of the so-called dwell tim o 18 ' 19 , 
g ^> 1 the normal-state conductance of the system, and d 
its mean level spacing. (However, the application of the 
method to other regimes requires only technical rather 
than conceptual modifications.) Here onwards we put 
h=l. 
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II. MAIN RESULT 

Deferring an outline of the technicalities of the anal- 
ysis to the last part of this work, we here merely an- 
ticipate that the correlator K(<f>i,<f>2) is obtained by a 
second-order perturbative expansion around the Usadel 
functional free energy. In this way we find, 



c r e 2 A 2 f°° 
Kr(<f>i,<h) = 5 — sin 0i sin 02 / dxi / dx 2 



y/xj + iy/xi + i 



sjx{ + cos 2 (0i/2) v'z 2 + cos 2 (0 2 /2) 
+ x /.x 2 + u/x 2 +cos 2 (0 1 /2))" 2 , (2) 



where 



cr = 



for highly transparent interfaces (T = 1) 
for weakly transparent interfaces (r <C 1), 

(3) 



and r is the transmission coefficient characterizing the 
transparency of the SN interfaces; see below. While 
Eqs. (|2|) and ((3|) are a new result for weakly transpar- 
ent interfaces, current fluctuations of the moderately 
weakly coupled quantum dot with highly transparent in- 
terfaces have been the subject of an earlier scattering 
theory analysis in Ref. [l2|. Within this approach both 
the diagonal contribution var J(0) = K((f), 0)^2- and the 
straightforward extension to the full correlation function 



K{4>\, 02) are represented as a double integral over eigen- 
values of the transmission matrix, which despite its dif- 
ference in form, is numerically identical. The agreement 
with these earlier results provides a check on the con- 
sistency of the field theory. Notice that the validity of 
the result is limited to values of the phase outside the 
domain 2 ^ |0 — tt\ < g^ 1 , a fact that was pointed out ear- 
lier in Ref. [12] • The main result of this work is that the 
Josephson-current fluctuations for the case of weakly and 
highly transparent interfaces differ by a universal factor 
of 2. This factor 2, in fact, can be anticipated by the 
observation that for small phase differences 0, one may 
establish contact to the theory of universal conductance 
fluctuations (UCFs) through a disordered quantum dot- 
i.e., for <C 1 



J(0) 



geA 



(4) 



where g is the (dimensionless) conductance through the 
disordered quantum dot. Using the fact that for the or- 
thogonal ensemble 2 ^ 



var g 



for high barriers (f < 1) 



| for transparent contacts (T = 1), 
one is led to the expectation (0 <C 1) 



(5) 



var J(0) 



16 

,2 A 2 



b 2 for high barriers (r -C 1) 
b 2 for transparent contacts (r = 1), 

(6) 
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in agreement with Eqs. |2]) and ([3]). 

Figure [1] shows the Josephson-current fluctuations 
^/var J((j)) = ^/A(</>, <j>) as a function of </>. We observe 
that the maximum falls to a fixed value (f> ~ 2.7. The 
strong increase of the fluctuations is a manifestation of 
the fact that <f> = it defines an instability of the Joseph- 
son action, in the vicinity of which minute changes of 
the disorder/boundary geometry trigger drastic changes 
in the configuration-specific J ((/))■ 



III. FORMALISM 

Having discussed our results for the SNS quantum dot, 
we next provide a brief outline of the formalism. 

Assuming M channels coupling the disordered quan- 
tum dot to the superconductors (with each M/2 modes 
propagating to the left and to the right, respectively^) 
the Hamiltonian of the SNS junction is given by 



H = Ha 



Hi 



H\ + H} + H 2 



where (i = 1,2) 

#d = £i^)ic;wvi, 



(7) 



(8) 



Hi = W dE\xi(E))H^{ X i,{E)\, (9) 

aa f 

H l = Hl dE(\^)Wrixi(E)\ + \xi(E))Wt a {t 



(j,a 



(10) 



{|Vv»)l denotes a set of dot states (for further conve- 
nience we use the Nambu formalism; i.e., the states 
comprise particle and hole degrees of freedom in a sin- 
gle object), and {\Xa(E))} a set of scattering states into 
superconductor i. The dot Hamiltonian is of the form 
^dot = Cg 11 (g) h where h is drawn from the orthogonal 
Gaussian ensemble with the variance set by the spectrum 
width A, £?B dG is the Boguliubov-de Gennes Hamiltonian 
of superconductor i with order parameters = Ae'* , 
and the coupling matrices W l are characterized by the 
transmission coefficients^ T, WiW* = v where 

f{T) = ^ — 1 — pVl — r and fi&li describes a diagonal 
matrix with entries 1 for the M/2 open channels con- 
necting the dot to superconductor i and otherwise (We 
assume that the sets of dot states coupling to supercon- 
ductors 1 and 2 are disjoint.) The capacitance of the 
dot is assumed to be sufficiently large that the Coulomb 
blockade can be ignored. 

To compute Josephson-current correlations we intro- 
duce a supersymmetric field integral representation of the 
free energy^: 



K 



{<h,<h) = (~) 2 / 



d(pidip 2 



I Vl=¥'2=0: 



(11) 



where 



with 



$) = J VQe 



-S[Q] 



(12) 



M / 
S[Q] = y Restrln (l 



/(T) 



Vuj 2 + A 2 



LOG 



ph 



Aph i$o 
&2 e 



Q 



(13) 



Here, Q is a 16-dimcnsional matrix acting on the prod- 
uct of particle-hole (ph) space, a two component space 
(f) required to distinguish between the two supercur- 
rents, a boson-fermion (bf) space implementing the su- 
persymmetric structure of the theory, and a time re- 
versal (tr) space required to correctly describe the be- 
havior of the system under time reversal. Further, 
the symbol "str" denotes the supersymmetric general- 
ization of the matrix trace, and lj = diag(cj 1,^2) and 
<J) = diag(</>i + ipi, (f>2 + ip%) are diagonal matrices in f- 
space where lo\^ and 4>i,2 oc l bf are the two energy and 
phase arguments entering the integral representation of 
the supercurrent and the sources tpi t 2 = diag(v?i i 2, 0)bf 
break supersymmetry. 

Referring for a detailed discussion of the derivation of 
Eqs. pip - (fT3"|) and of the internal structure of the ma- 
trix Q to Ref . fl(| . we here merely recapitulate, that it 
originates from the coherent state representation of the 
free energy along the standard procedure, including the 
generalization of the partition function to a supersym- 
metry formulation, RMT ensemble average, Hubbard- 
Stratonovich transformation and expansion in Goldstone 
modes around the metallic saddle point <7g h . 

Notably, the nonlinear constraint Q 2 — 1 indicates 
that the matrix Q represents a generalization of the 
Green function g central to the quasiclassical approach. 
Indeed, it is straightforward to see that a variation of Eq. 
(fT3"| (under the restriction Q 2 = 1) leads to 



, w(T? h + A cos $er.? h 



0, 



(14) 



i.e., an equation which, upon identification of the ma- 
trix Q with the quasiclassical Green functions, imme- 
diately is recognized as the zero-dimensional limit of a 
Usadel equation (extended, however, to a larger struc- 
ture accommodating more than one impurity averaged 
observable.) Void of elements coupling between the two 
observables (f-space), Eq. (fl4|) admits a block-diagonal 



solution Q — diag(gi, 32) ® 1 



bf = RafR- 



where 



fji 



(oJiof + A cos $ lCT P h ) (15) 

-I- A 2 v / 



are the quasiclassical Green functions computed for phase 
difference $1.2 and the second representation expresses 
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the solution as a rotation away from the metallic ref- 
erence point, (Tg 11 . The ensemble-averaged Josephson- 
current may now be calculated from the saddle-point ac- 
tion 



where we introduced Dk(4) := ^gggp i . Perform- 
ing the twofold derivative with respect to the sources 
followed by Gaussian integrals over the W's one obtains 
the final result, Eqs. © and ©. 



■S[Q] = -^-strln (l + / 2 (r) + 2/(r) 



V^ 2 + A 2 cos 2 $ 



Vw 2 + A 2 



(16) 



and one obtains for highly transparent interfaces 

K£ roc 

J{4>) — -7, — A sin (f> / dx 
2tt Jq 

1 1 

^/x 2 + 1 + ^Jx 2 + cos 2 <f>/2 ^Jx 2 + cos 2 0/2 ' 

(17) 

while for weakly transparent interfaces 



rri . eMI\ . 
J(<p) —— — Asm< 



dx- 



COS* 



0/2 



1 

(18) 



Equations (fTTll and (fT8|) coincide with quasiclassical 
resultsli^^^. 

Josephs on- current correlations can now be explored by 
defining*^ 



(19) 



where the generalized rotation matrix T describes fluctu- 
ations around the Usadel mean field and the free energy 

= JVTeM-SlQ})- 
Going beyond the mean-field level by substituting the 
generalized parametrization, Eq. (|19p . into the field in- 
tegral, one finds that fluctuations around T = 1 are pe- 
nalized by a parameter cos(<pi/2)g which is much larger 
than unity, unless \<f> — tt| < g _1 lies in the anomalous 
window discussed above. While in the latter case the 
correlator, Eq. (fTTj) . has to be evaluated by full integra- 
tion over the manifold of T's (cf. the analogous situation 
with computing the fine structure of spectral correlation 
functions on energy scales smaller than the single parti- 
cle level spacing^), in general it is sufficient to pertur- 
batively expand the matrix T = t + W + \ W 2 + . . . to 
second order around unity. Expansion of the logarithm 
then leads to the Gaussian actions for highly and weakly 
transparent interfaces, 



TM 



4 2 2iW = Y str 



4 

M 



(20) 



1 , ; , A ~ sin *o% h ofW 

1 + Vlo 2 + A 2 2 3 



(21) 



IV. GENERALIZATION 

While the outline above was specific to the case of 
a spatially structureless quantum dot, the generaliza- 
tion to more complex geometries is straightforward. Un- 
der conditions where the spatial variation of the qua- 
siclassical Green function is no longer negligible, Eq. 
(fT3| generalizes to the action of a matrix field Q. For 
example, for an extended diffusive N region (i.e., for 
E c < A, Eg) and weakly transparent interfaces, S[Q] = 
-(nDu/8) J d d rsti(dQdQ) + S4Q} + S c [Q], assumes the 
form of a generalized diffusive a model,— where S u [Q] = 
— (■kv/2) J d d r str(Q6jo-^ h ) and S C [Q) describes the cou- 
pling to the superconductor. Similarly, for a (nearly) 
clean system, S[Q] = -v F J d d rdnstr(To-fn ■ dT~ v ) + 
Su[Q] + S C [Q], becomes the free energy of the ballis- 
tic er-modcl 25 . In complete analogy to our discussion 
above, the variation of these actions leads to the gen- 
eral Usadel^ and Eilenberger— equations, respectively. 
To compute Josephson-current correlations, one again 
employs the parametrization, Eq. (|19[) . where, however, 
both R and T are field configurations with nonvanishing 
spatial variation. By evaluating the action of a general- 
ized quadratic variation T = 1 + W + \ W 2 and subse- 
quent computation of the correlator, Eq. pip , one can 
then, in principle, obtain the Josephson-current corre- 
lations of any SN system category (if amenable to the 
approximation schemes of quasiclassics.) However, in 
cases where the solution of the Usadel and Eilenberger 
equations displays a complex spatial profile, the concrete 
computation of the generalized Gaussian integrals over 
W can be cumbersome. Under these conditions one ex- 
pects the result for the current correlations to be less 
universal (e.g., dependent on the system geometry, dis- 
order concentration, etc.) than in the quantum dot case 
discussed above. 



V. SUMMARY 

In summary, we showed that the mesoscopic 
Josephson-current fluctuations through a weakly coupled 
disordered quantum dot differ by an universal factor of 
2 for the cases of weakly and highly transparent inter- 
faces. In our calculations we applied a supersymmet- 
ric field-theoretical approach. Although we concentrate 
on the case of short junctions (and use RMT to model 
the normal metal), this method allows for a microscopic 
derivation capable to cover a wide range of parameters 
characterizing the SNS junction. The observation of the 
Josephson-current correlations should be feasible using 
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